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Abstract. In this paper the Boltzmann equation describing the carrier transport 
in a semiconductor is considered. A modified Chapman- Enskog method is used, in 
order to find approximate solutions in the weakly non-homogeneous case. These 
solutions allow to calculate the mobility and diffusion coefficients as function of 
the electric field. The integral-differential equations derived by the above method 
are numerically solved by means of a combination of spherical harmonics functions 
and finite-difference operators. The Kane model for the electron band structure is 
assumed; the parabolic band approximation is obtained as a particular case. The 
numerical values for mobility and diffusivity in a silicon device are compared with 
the experimental data. The Einstein relation is also shown. 
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1. Introduction 

Many commercial simulators of microelectronic devices are based on the well-known 
Drift-Diffusion equations. In this model the electron current density J n is given by 
the equation p|, H 

(1) J n = fi n pE + A»V x p, 

where p is the electron density, E the electric field, p n and D n the mobility and the 
diffusion coefficient, or diffusivity, respectively. The symbol V x denotes the gradient 
operator with respect to x. For low electric field the two transport coefficients p n 
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and D n are related by the well-known Einstein relation 
(2) D n = (M^) „ n 

where ks is the Boltzmann constant, Tl the lattice temperature and e the absolute 
value of the electron charge. The transport coefficients fi n and D n are often assumed 
constant. This approximation becomes inadequate in presence of high electric fields. 
In these cases Eq. ([I]) remains valid only if the transport coefficients are considered 
as functions of the electric field. Many expression for these functions were proposed 
(see, for instance, Refs. |§, f|, They were obtained by fitting experimental 

data or Monte Carlo simulations [|J, [12|]. In this paper we obtain the functions 



/j> n {E) and D n (E) directly from the Boltzmann transport equation. 
The paper is organised as follows. In Sec. 2 we briefly introduce the Boltzmann 
transport equation for an electron gas in a semiconductor. Only electron-phonon 
scatterings are considered. So that, electron-electron and electron-impurity interac- 
tions are assumed negligible. In Sec. 3 we perform a Chapman-Enskog expansion, 
and in Sec. 4 a numerical scheme is proposed in order to solve the obtained equa- 
tions. In the last section we show numerical results for fi n and D n , which we compare 
with experimental data. Further, we analyze the validity of the Einstein relation for 
different values of the electric field. 

2. Basic equations 

For an electron gas in a semiconductor device the coupled system Boltzmann-Poisson 
equations f2[ writes 

(3) ^ + v(k) ■ V X F - ■ V k F = Q(F), 

(4) V x ■ E = Z - (n d {x) - N A (x) -jFdk^j 

where the unknown function F(t, x, k) represents the probability of finding an elec- 
tron at the position x, with wave-vector k, at time t. The wave-vector k belongs 



to M 3 . The integrals with respect to k are performed over the whole space and 
the parameter h is the Planck constant divided by 2tt. The group velocity v(k) 
depends on the band structure and it will be defined in the following. The symbol 
Vk denotes the gradient with respect to the variables k. The functions Na , -/Yd are 
the concentration of acceptors and donors and e is the dielectric constant. Here we 
assume that the low-density approximation holds, so that Q is a linear operator. If 
K(k, K) is the symmetric scattering kernel, then 

(5) Q(F)(t, x, k) = J K(k, K)F(t, x, k*) dK - F(t, x, k) J K(K, k) dK ■ 

We consider the scattering kernels describing acoustic phonon interactions (in the 
elastic approximation) 

and non-polar optical phonon interactions 

K op (K, k) = (^ n + l \g X5(e(K) - e(k) ± tiw) . 

Here, Xq and X are constant and the overlap factor Q for the conduction band is 
taken equal to 1 (see J4|, flO||). We consider negligible the ionized impurity scattering; 



i.e. we assume low density of doping. With these assumptions, the collision operator 
esplicitelly writes 

Q(/)(t, x, k) = (n, + 1) / X5(e(K) - e(k) - hu)f{t, x, k») dK 
+ rig J X6(e(k*) - e(k) + hcu)f(t, x, k,) rfk* 

(6) + J X S(e(K) - e(k))/(t, x, k.) dK - ?(k)/(t, x, k) , 
where 

P(k) = n g / X5(£(k*) - e(k) - M dk, + (n g + 1) 

(7) x J X5(e(K) - e(k) + ^) dK + / X 5(e(k*) - e(k)) dK ■ 



For every admissible function F(t, x, k), we have [7j 

J Q(F) dk = 0, 

which implies mass conservation. The band structure is assumed to be spherically 
symmetric. Hence, the electron energy e is given by the equation 



(8) ^> = 2^' 

where m* is the reduced electron mass. We choose the Kane formula 

(9) 7 (e) = e(l + ae), 

which is simple but sufficiently accurate to describe high field regime. By putting 
a = one obtain the usual parabolic band approximation. The velocity v(k) is 
explicity given by 

1 frk 

(10) v(k) := -V k e(k) 



h m*{2ae + 1) 



3. The Chapman-Enskog expansion 

The Chapman-Enskog method || allows us to find approximate solutions of the 
Boltzmann equation. They are valid only for a very small space domains. More- 
over, these regions must be far enough from boundary layers. Nevertheless, these 
solutions have a great importance since they allow to obtain the transport coeffi- 
cients. In the framework of electron transport in semiconductors, we expect that 
such solutions are valid in spatial domain, where the electric field and the doping 
are almost constant. Therefore, they should not furnish good results near junctions 
or boundaries. 

We follow the standard scheme p] , assuming that the solution of Eq. (|3]) is approx- 
imately expressible by the following relation 

(11) F(t, x, k) ~/(t,x,k) +<*/(*, x,k). 
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As usually, the typical Chapman-Enskog constrain 

(12) J 5f(t,x,k)dk = 

is imposed. Therefore using Eqs. (PT])-(0) we obtain that the electron density is 

p(t,x) := J F(£,x,k)dk~ J /(*,x,k) dk. 

We assume that time and space partial derivatives of f are of the same order of 5f 
and that the corresponding derivatives of 5f are negligible. By inserting Eq. (|ll|) 
into Eq. @ and splitting terms of different orders, we obtain the following equations 

(13) ~E.Vk/ = Q(/), 

df c 

(14) -± + v(k) • V x / - • V k 5/ = Q(5f) . 

In this approximation, the Poisson equation becomes 

(15) V x • E ~ - e (n d - N A - J fdkj . 

In order to solve Eqs. (pf)- fljjj) we assume that E is constant. This is a reasonable 

assumption, if the difference between the density of electrons and doping is negligible 

in the small domains, where we consider the Boltzmann equation. 

Of course, we solve before Eq. ([T3| ) and then the next one. Since Eq. ( |I"3l) is linear 

and it does not depend explicitly on the variables (t, x), a solution can be written 

as 

(16) /(t,x,k)=p(t,x)s(k). 

We note that the function g(k) is also the stationary homogeneous solution of Eq. @ 
for constant electric field, verifying the condition 

(17) Jg(k)dk=l. 

Now, we consider Eq. (|TJj). Using flrp|), it becomes 

dp 



9t+ v(k).VxP 



^E-V k 8f = Q(8f). 



By integrating Eq. (|T^) with respect to k, we find a compatibility condition, as 
usually arises in integral equations. Assuming the reasonable hypotesis / V^Sfdk = 
0, we obtain 

(19) ^ + V -V x p = 0, 
where 

(20) V := J g{k)v{k)dk 

is the constant macroscopic electron velocity in the stationary homogeneous case. 
This vector is parallel to E (see appendix A). By eliminating the time derivative of 
p in Eq. (|i"8|) by using Eq. (|19"D, we get 



(21) g(k) [v(k) - V] ■ V xP - |E • V k 8f = Q(5f) . 
The form of Eq. (^T|) suggest us to assume 

(22) 5/(t,x,k) = |V x p(t,x)|/i(k). 



Now, Eq. ([12] ) implies that 

(23) Jh(k)dk = 0. 

We denote by u the unit vector in the direction of V x p; i.e. V x p = |V x p|u. It in 
general may depend on the variables (t, x). Using Eq. (^), it is easy to see that 
Eq. (^l|) becomes 

(24) g(k) [v(k) • u - V • u] - |E • V k /i = Q(h) . 

n 

The case u parallel to E is the most meaningfull. 
4. Approximate equations 

In order to obtain / and Sf we have to solve the following set of equations for the 
unknowns g and h 

(25) -|E-V k <7 = Q(<?), 



(26) <7(k) [| v(k) | cos 9 - V] - |E ■ V k /i = Q(h), 

(27) Jg{k)dk=l , Jh{k)dk = 0. 

Here # is the angle between v(k) and u, and V = |V|. 

Since g and fa are functions of the 3-dimensional variable k, we use a spherical har- 
monics expansion to reduce the dimension of the space of the indipendent variables. 
Taking in account the symmetries of the problem and considering only the first two 
terms of the expansion, we can choose 

(28) g(k)~go(e)+ gi (e)cos9, 

(29) h{k) ~h {e) + h 1 (e)cos6. 

The use of the Galerkin method allows to derive from Eqs. (^5|)- (|26|) a set of equa- 
tions for go,gi,h and hi. This approach recalls the well-known method of the 
spherical harmonics expansion to solve the Boltzmann-Poisson system, ||, ||, 
Since Eqs. contain the same diffusion and collisions terms of the Boltz- 

mann Equation, many calculations are performed following Ref. ||. Then, the 
system (p5|)-(p6|) gives the following set of ordinary differential-difference equations 



(30) - t E(^f±g x + g>\^- = Q x {g ) 

\j(e) ) v{e) 

(31) -e#<7o4r = <3 2 (<7i) 

v(e) 

(32) -Vgo(e)^ + 9l (e) - zE (-^hi + h[) JL = Qi{h ) 

v{e) \j(e) J v{e) 

(33) -Vgi{e) -L + g {e) - zEh' = — -r— = Q 2 (h) , 

v(e) v{e) 

where the prime indicates the derivative with respect to e, 
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and 



Qi((p) = (n 9 + 1)X a(e + Tuo)ip{e + hu) + n q X o(e — fuu)(p(e — hu) 

(34) - [n q X a{e + hu) + (n, + 1)X a{e - <p(e) , 

(35) Q 2 (y?) = - [u q X a(e + hu) + (n g + 1)3C <x(e - Uu) + X <r(e)] (p(e) . 



The density of states cr(tt) is given by 

: = y 5( e (k) - u) dk = AV2n ^- 
where H(u) is the Heaviside step function. 



'm 



H(u)^/j(u)^(u), 



5. Physical model and numerical results 

We use the following values for the parameters, appropriate for a silicon device: 



m* = 0.32 m e 


T L = 


300 K 


= 0.063 eV 


x - {DtK)2 

87T 2 pUJ 


D t K 


= 11.4 eV A" 1 


p = 2330 Kg m~ 3 


k B T L „ 2 


£<f = 


9 eV 


vq = 9040 m sec -1 . 


a = 0.5 eV- 1 









In the table, m e is the electron rest mass. We have solved numerically Eqs. (p0|)- 
(33) with the conditions for different values of the electric field. The numerical 
procedures are similar to that used in Ref. ||. The kinetic definition of the current 
density is 

(36) J n := J F(t,x,k)v(k)dk~ p(t,x)V +\V x p(t,x)\J h(k)v(k)dk. 



Since V = Vu and 



/i(k)v(k) dk ~ / [h (e) + h x (e) cos#]v(k) dk 



(37) 



^i( £ ) V(k) i ; U v(k) dk = u / W V(k) / U)2 dk . 



v k 



V k 



a comparison between Eqs. (p6|)-(p^) and (JT]) gives 

(38) ^„ = ^ and D n = ( h x {e) (v f k) ' dk. 

|E| J | v (k)| 

They are functions of E, because V and hi depend on this parameter. 

In fig. 1 we show the mobility /i n against the electric field. The values are compared 

with two curves which fit experimental data by means of simply formulas (see [|ll[ 



pp. 94-98). For low values of the electric field there is a small difference. It is 
meaningless because many different values (up to 20 % ) of low field mobility are 
reported in literature (see Jill, PP- 81-82). In fig. 2 the value of diffusivity is shown 
together with experimental data (see p. 6715). The agreement is very good for 
all range of |E|. Fig. 3 and fig. 4 compare the mobility and the diffusivity using Kane 
model and the parabolic band approximation. The last figure shows that Einstein 
relation is valid only for moderate fields. 
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Appendix A 

The unknown g must satisfy Eq. (|25|) , which contains only two vectors, k and E. 
Since g is a scalar function, the only possibility is that g depends on k only throught 
the scalars |k| and k ■ E i.e. g = g(|k|, k ■ E). Therefore, from Eq. (|20|) we have 

V:=/ 9 (k)v(k) (flt = / 5 (|k|,k.E)-^ TTjd k. 

Hence, it is evident that V is parallel to E. 
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Figure 1: Mobility as a function of electric field: dots indicates values obtained using 
Kane band model, solid line a fit with parameters of Canali et. al. and dashed line another 
fit with parameters of Caughey and Thomas (see Ref. |Tl|| ). 
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Figure 2: Diffusivity as a function of electric field: solid line indicates values calculated 
using Kane band model and squares and circles experimental data (see Ref. pi]). 
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Figure 4: Diffusivity versus electric field: Kane band model (solid line) and parabolic 
band approximation (dashed line). 
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Figure 5: Comparison between diffusivity calculated directly (dashed line) and using 
Einstein relation (solid line), both with Kane band model. 
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